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ABSTRACT 

We make a more general determination of the inflationary observables in the 
standard 4-D and 5-D single-field inflationary scenarios, by the exact reconstruc- 
tion of the dynamics of the inflation potential during the observable inflation with 
minimal number of assumptions: the computation does not assume the slow-roll 
approximation and is valid in all regimes if the field is monotonically rolling down 
its potential. 

We address higher-order effects in the standard and braneworld single-field 
inflation scenarios by fitting the Hubbble expansion rate and subsequently the 
inflationary potential directly to WMAP5+SN+BAO and Planck-like simulated 
datasets. 

Making use of the Hamilton- J acobi formalism developed for the 5-D single-field 
inflation model, we compute the scale dependence of the amplitudes of the scalar 
and tensor perturbations by integrating the exact mode equation. The solutions 
in 4-D and 5-D inflation scenarios differ through the dynamics of the background 
scalar field and the number of e-folds assumed to be compatible with the obser- 
vational window of inflation. 

We analyze the implications of the theoretical uncertainty in the determina- 
tion of the reheating temperature after inflation on the observable predictions 
of inflation and evaluate its impact on the degeneracy of the standard inflation 
consistency relation. 

We find that the detection of tensor perturbations and the theoretical un- 
certainties in the inflationary observable represent a significant challenge for the 
future Planck CMB measurements: distinguishing between the observational 
signatures of the standard and braneworld single-field inflation scenarios. 



This work have been done in the frame of Planck Core Team activities. 
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Subject headings: cosmology: cosmic microwave background, cosmological pa- 
rameters, early universe, observations 



Introduction 



The primary goal of particle cosmology is to obtain a concordant description of the 
origin and early evolution of the Universe, consistent with both unified field theory and as- 
trophysical and cosmological measurements. 

Inflation is the most simple and robust theory able to explain the astrophysical and cos- 
molog ical observations, p roviding at the same t i me self-consistent primordial initial condi- 



Linde 1983 


) and the mechanisms for c 


quantum generation of the scalar and tensor per- 


turbations ( 


Mukhanov & Chibisov 1981: 


Hawkine 1982: 


Starobinskv 1982: 


Guth & Pi 1982: 



Bardeen et al. 19831 : lAbbott fc Wise 1984( 1. 

In the simplest class of inflationary models, inflation is driven by a single scalar field 
(inflaton) with some potential V(<f)), minimally coupled to Einstein gravity. The per- 
turbations are predicted to be adiabatic, nearly scale-invariant and Gaussian distributed, 
resulting in an effectively flat Universe. At the leading-order in s low-roll approximation 
( ISteinhardt fc Turner 1984 : ISalopek fc Bond 19901 : iLiddle et al. 1994 ) the amplitudes of scalar 
and tensor perturbations on a specified comoving wavenumber k, are related through the 
consistency equation: 

~2~ 



A 5 



where: A? s oc k ns ~ l and A\ oc k nT are the amplitudes of scalar and tensor perturbations 
respectively and rig and are their tilts. The consistency equation may be regarded as 
an independent test of single-field inflationary models as it does not depend on the specific 
functional form of the inflation potential. 



Recent WMAP 5-year CMB mea surements, alone (IDunklev et al. 20091 ) or complemented 
with other cosmological datasets ( Komatsu et al. 20091 ). support the standard inflationary 
predictions of a nearly flat Universe with adiabatic initial density perturbations. 
In particular, the detected a nti-correlations be tween temperature and E-mode polarization 
anisotropy on degree scales (INolta et al. 20091 ) provide strong evidence for correlation on 
length scales beyond the Hubble radius. 

Despite the successes of inflationary cosmology, recent p roposals in theoretical p hysics 
motivated by the developments in superstring and M-theory (IHorava fc Witten 19961 ). sug- 
gest that our four- dimensional Universe could lie on a brane embedded in higher-dimensional 
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space-time (see e.g. Rubakov 2001, Maartens 2004 and referenc es therein). 



In particular, in the type II Randall- Sundrum model (RSII) ([Randall fc Sundrum 1999a] ; 
Randall &: Sundrum 1999bl ) our four- dimensional (4-D) Universe is a brane with positive 
tension A embed in a five-dimensional (5-D) anti-de Sitter space-time (AdSs). At sufficiently 
low energies (p << A) the standard cosmic behavior is recovered and the primordial nucle- 
osynthesis constraint is satisfied, provided that A >(1 MeV) 4 . 

The simplest way to realize inflation in R SII model is to have a single scalar field confined 
to the brane and only gravity in the bulk (IMaartens et al. 20001 ) . In this case the Fiedmann 
equation is modified at so that Hubble parameter H oc p rather than H oc p 1 / 2 as in 4-D 
case, leading to significantl y modifications of the amplitudes and scale dependencies of scalar 
and tensor perturbations (IBinetruy et al. 20001 ). 



The observational constraints on inflationary parameters in 5-D scenario, made in gen- 
eral by using the slow-roll (SR) approximation in the high energy regime (p >> A), show 
that to the leading order in SR the consistency equation has precisely the same form as in the 
standard 4-D scenario , the relationship between inflationary observab les being independent 
on the brane tension (iTsujikawa fc Liddle 20041 : ISeery fc Taylor 20051 ) . 
The degeneracy of consistency equation was associated with the fact that 5-D inflationary 
observables smoothly approach their 4-D counterparts as the brane decouples from the bulk 
approaching the low energy regime (p << A). The main assumption made by these works 
is that the back-reaction due to metric perturbation in the bulk can be neglected. This 
assumption is valid to the leading-order in slow-roll approximation, as the coupling between 
inflation fluctuations and metric perturbation vanishes. 



Recently it was shown (IKoyama et al. 2004 : iKoyama et al. 2005at iKoyama et al. 2005bl ) 
that the sub-horizon inflation fluctuations on the brane excite an infinite ladder of Klauza- 
Klein modes of the bulk metric perturbations to second-order in slow-roll parameters. If the 
back-reaction is take into account, the amplitude of the scalar pertur bations receives second - 
order slow-roll corrections in ad dition to Stewart-Ly th corrections (IStewart fc Lyth 19931 ). 
of the same order of magnitude (IKoyama et al. 20081 ). The degeneracy of consistency equa- 
tions doe s not hold when the second-order corrections in SR expansion for perturbation s are 
included rtCalcagni 20031 : ICalcaeni 20o3 iRamirez fc Liddle 2004! : ISeerv fc Tavlor 2005^ . 



One of the most anticipated results of forthcoming high precision CMB experiments is 
probing the physics of inflation and in particular the reconstruction of the inflation potential. 
In order to have a robust interpretation of upcoming observations it is imperative to under- 
stand how the reconstruction process may be affected by the degeneracy of the inflationary 
observables. In this paper we aim to make a more general determination of the inflationary 
observables in 4-D and 5-D inflationary scenarios, by exact reconstruction of the dynamics of 
the inflation potential during the observable inflation with minimal number of assumptions 
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( jLesgourgues at al. 20081 ; lHamann et al. 20081 ) 



?aking the advantage of the formalism c 


eveloped for the standard single-field inflation 


Peiris et al. 2003 




Peiris & Easther 2006a 




Peiris & Easther 2006b 




Martin & Rineeval 2006 



Lesgourgues at al. 20081 ; lAlabidi fc Lidsey 20081 ). we carry out similar calculations for 5-D 
inflation models by fitting the Hubbble functio n, H(<t>), and subsequently the inflationary 
potential, V (</>), directly to WMAP 5-year data JPunklev et al. 20091 : Ecomatsu et al. 2009h 
complemented with geometric probes from the Type la supernovae (SN) distance-redshift 
relation and the baryon acoustic oscillations (BAO) measurements and Planck-like CMB 
anisotropy simulated data. 

Our specific goal is to address higher-order effects in the standard and braneworld single-field 
inflation models and to analyze the sensitivity of the present and future CMB temperature 
and polarization measurements to discriminate between them. 

The paper is organized as follows. In Section 2 we review the Hamilton- Jacobi for- 
malisms for 4-D and 5-D sigle-field inflation models. In Section 3 we compute the scalar and 
tensor perturbation spectra for standard and braneworld single-field infation models by using 
the exact mode equation. In Section 4 we present the implementation of the Markov Chain 
Monte-Carlo methodology and describe the datasets involved in our analysis. Section 5 is 
dedicated to the analysis and the interpretation of our results: we present the derived bounds 
on the inflationary parameters, Hubble Slow-Roll parameters and the magnitude, slope and 
curvature of the infationary potentials obtained from the fits of 4-D and 5-D single-field 
inflation models to our datasets and analyze the possibility do disentangle between standard 
and braneworld scenarios by using the future Planck high precision CMB measurements. 
In Section 6 we draw our conclusions. 

Throughout the paper m 4 and m 5 denote the corresponding 4-D and 5-D Planck mass 
scales and we have set Gm\ = h = c = 1. Also, we denote by dot the derivative with respect 
to the time and by prime the derivative with respect to the scalar field. 



2. The four- and five-dimensional Hamilton- Jacobi formalism 

2.1. The 4-D single-field inflation case 

The Hubble Slow-Roll (HSR) formalism for the standard single-field inflation was set 
down in detail by Liddle, Parsons & Barrow (1994). 
The Friedmann equation in zero-curvature Universe is given by: 
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where: H = a/ a is the Hubble parameter, a is the cosmological scale factor, p = V + 
2 /2 is the total energy density, where V(<f)) and 4> 2 /2 are the potential and kinetic energy 
density terms respectively. Since the da rk energy contribution is strongly suppresse d by the 
exponential expansion during inflation (IMaartens et al. 20001 ; lLanglois et al. 20011 ) . we set 
to zero the dark energy term in the above equation. 
The equation of motion for the scalar field is given by: 



+ 3H<j) = -V . 



(3) 



Eqs. ([2]) and ([3]) can be written in the Hamilton- J acobi form, allowing to consider inflation in 



terms of H(<j)) ra t her than V(6) (ILiddle et al. 1994 : iKinnev 2002 : Easther & Kinnev 2003 
Peiris et al. 20031 : iKinnev et al. 20041 ): 



_±Ltf(0a(0) 



rn 



12tt 



ml 



H\<f>) 



4 

47T 

32tt 2 



m 



x v{4>) . 



(4) 

(5) 
(6) 

and (151) allow to 



For any value of H{(f>) Eq. j6]) can be used to find V(4>) while Eqs. 
convert 0-dependence into time-dependence. 

In the standard 4-D inflation the first three HSR parameters are given (ILiddle et al. 19941 ): 

ml H' 2 (4>) 



Vh{<P) 

CM 



4tt H 2 {<f)) ' 
ml H"{(j>) 
4^ H{(j>) ' 
m\ H\(j))H'"((P) 



(7) 



(9) 



16vr 2 H((f>) 

The dependence of V(4>) on H((f>) can be obtain by substituting en into Eq. (jHJ) leading to: 

1 



in 



3ml 



V{4>) = h\4>) 



3^(0) 



(10) 



The HSR formalism ensures that the condition for inflation to occur is precisely en < 1 and 
inflation ends exactly when en = 1. 
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2.2. The 5-D single-field inflation case 

In the 5-D inflation case the Eq.(J2]) receives an additional term quadratic in energy 
density: 



where: p = V + <f)/2 is the total energy density and A is the brane tension. The scalar field 
cf) is assumed to obey the same equation of motion as in 4-D standard inflation as given by 
Eq-©. In the low-energy regime (p << A) the quadratic term in Eq. ffTTl) can be neglected 
and one recover the behavior of the 4-D standard cosmology. In high-energy regime (p >> A) 
the deviation from the standard expansion changes the amplitudes and scale-dependence of 
cosmological perturbations. 

Hereafter we will make use of the approach developed by Hawkins and Lidsey who derived 
a general formalism for 5-D inflation case valid in all regimes, h aving many of the prop 



erties of the Hamilton- J acobi formalism in 4-D standard inflation ( Hawkins fc Lidsey 2001 



Hawkins fc Lidsey 20031 ). They defined a quantity y(<p) with the same role as H(<p) in the 



case of 4-D standard inflation: 



with the inverse relation given by: 



2Ay 2 (0) 

P = ~\ WT\- 13 

1 - y 2 ((j)) 



In terms of y(4>) the Friedmann equation ( fill reads as: 



r2 16ttA y 2 (cj>) 



3m| (1 - y A {(p)Y 

where the restriction y 2 < 1 is imposed, implying that y(<p) is proportional to H(<f>) in the 
low-energy limit, y — > (p/A — > 0). 

The Ham ilton- Jacobi equations, analogues to Eqs. (|4D - (j6"D for 4-D standard inflation are 



given by (IHawkins fc Lidsey 20031 : iRamirez fc Liddle 20041 ) : 



47T 

y'(0)a'(0) = 2^)a(0) , (15) 

mi 

\m\ x 1 



, i , , ( 16 ) 
Sir J i — y 

tju a\ 4tt (1 + y 2 ) • 

#(0) = 2 71 2T^' ( 17 ) 

m| (1 - y 2 ) 
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and the dependence of V(4>) on y{<p) can be obtained by combining Eqs. (jl6p and ffTTj) leading 
to: 

2Ay 2 



7(0 



y 



6tt \1 — y z 



The first three HSR parameters in terms of y(4>) reads as (IRamirez &: Liddle 2004 ): 



?7h 



Am 2 
A77I4 

H 2 c 



1/2 



1/2 



/2 



1 + ?/ 2 ) # 3 ' 

y if" 



% 3 



[l+y 2 )H 2 (l + y 2 ) 3 H 3 



Vh 



(18) 

(19) 
(20) 
(21) 



The above definitions of HSR parameters are valid in all regimes, generalizing the previous 
ones, preserving at the same time many of the inflation key properties: they are obtained by 
demanding the condition for inflation to occur precisely for e# < 1 and to end exactly when 
6h — 1. Also, they are preserving the lowest-order slow- roll definitions of the scalar spectral 
index, ns — 1 — 4e# + 2t}h, and of its running, dns/dlnk = 5enVn — 4e^ — 2^. 



3. The four- and five-dimensional exact mode equation 

The scale dependence of the amplitudes of the scala r (S) and tensor (T) perturbation s 



can be exactly obtained by integrating the mode equation (jMukhanov 19851 : iMukhanov 19891 ) 



4+ (k 2 -^)u k = 0, (22) 

V Z S,T/ 

where primes denote the second derivatives with respect to the conformal time. 
The numerical evaluation of the spectra involves solving Eq. (l22]) for each value of the 
wavenumber k, the evolution of \uk\/zs,T to a constant value defining the observable power 
spectra Vs,t- The solutions differ through the evolution of the background scalar field and 
the prior on the number of e-folds assumed to be compatible with the observational window 
of inflation. 



3.1. The 4-D single-field inflation case 



We compute the amplitudes of scalar and tensor perturbations by using the standard 
inflation numerical module from Lesgourgues et al. (2008). For each wavenumber k in a given 
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range the code integrates Eq. (1221) in an observational inflationary window corresponding to 
a number of AN e-folds, imposing that k grows monotonically to the wavenumber k* that 
leaves the Hubble radius when = 0*, eliminating at the same time the models violating 
the condition for inflation (e# < 1). 

For the purpose of present analysis we reconstruct the Hubble expansion rate H((f> — 0*) 
from the data by using the Taylor expansion up to the cubic term: 

H{4> - 0,) =H* + H:(<P- 0,) + \H'l{<j> - 0*) 2 + - 0*) 3 , (23) 

2 o 

equivalent to keeping the first three HSR parameters. We consider wavenumbers in the range 
[5 x 10 -6 — 5] Mpc -1 needed to numerically derive the CMB angular power spectra and the 
Hubble crossing scale fc* = O.OlMpc -1 . 

The analysis however depends on the prior on the interval over which the dynamics of the 
background field is tracked. Actually, the standard inflation numerical module integrates 
Eq. (1221) from the time at which k/aH = 50 until dlnPs,T/ 'dXna < 3 x 10~ 3 . This choice 
ensures that inflation started enough time before the observational range and ends (en = 1) 
enough time after the smallest observable scale leaves the Hubble horizon scale /c* = 0.01 
Mpc -1 , leading at the same time to an accuracy of ~ 0.1% i n final power spectra amplitudes, 



that is smaller than the expected sensitivity of CMB data (IHamann et al. 20081 ). For these 



reasons we choose to keep this time integration window for our comp utation 



The power spectra of scalar and tensor perturbations are obtained as (jStewart fc Lyth 1993 
Copeland et al. 19941 ): 



^ .,. k 3 ( H\ \u k \ 2 ^ ... 16/c 3 \u k \ 2 

2n z \ J a z ixm\ a 2 

where z$ = a<p/H for scalars, = a for tensors and the temporal evolution of the scalar 
field is given by Eq.(j2]). 



3.2. The 5-D single-field inflation case 

There is supporting evidence for the use of the exact mode Eq. (l22l in braneworld con 



text, because its derivation does not involve the Friedmann equation (ILiddle fc Lyth 20001 ) . 
We modify the standard inflation module from Lesgourgues et al. (2008) to compute the 
power spectra of scalar and tensor perturbations for the single-field braneworld inflation. 
As in the case of 4-D standard inflation, the Hubble expansion rate H(<fi — 0*) is obtained 
from the data by the Taylor expansion up to the cubic term the neighborhood of the pivot 
scale K=0.01 Mpc -1 . 
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For each wavenumbers k in the range [5 x — 5] Mpc -1 we integrate Eq. fl22|) keeping the 
same time integration window as in the 4-D inflation case. In this way, the identical scales 
(wavenumbers) encompass the same number of e-folds at Hubble radius crossing fc* in both 
4-D and 5-D inflation ary scenarios ensuring the same accuracy in the reconstruction of the 
inflationary potential (IKinney 2002l ). As in the 4-D case we impose the condition that each 
mode k grows monotonically to the wavenumber k* and we eliminate those models violating 
the condition for inflation (e# < 1). 

Taking zs = aH/<fi for scalars, z? = a for tensors and the evolution of the scalar field 
given by Eq.(fT6l). the power spectra of scalar and tensor perturbations are then obtained as 
jRamirez & Liddle 2004 : iKovama et al. 2008h : 

2 



V s (k) 



I: 



25vr 5 



H 



[1 + /C(/3 2 )] , V T (k) 



16k 6 \u k \ 



(25) 



The correction JC(P 2 ) to V s(k) is solely due to the coupling between the inflation field and the 
bulk metric perturbations (IKovama et al. 2005al ; IKovama et al. 2005bl; IKovama et al. 20 08). 
For each wavenumb er k we obtained lC(/3 2 ) by numerical computation ( Koyama et al. 20081 ). 
taking (3 2 given by (IKovama et al. 20081 ): 



P 2 



1+ T7 



h) 



-1/2 



(26) 



One should note that, although the power spectra of the tensor perturbations in 4-D and 
5-D inflation have the same form, they differ through their dependencies on the cosmological 
scale factor a(0) and on the conformal time: rj(a) = J da/a 2 H. 

Defining the amplitudes of scalar and tensor power spectra as ( Cop eland et al. 1994 )^1: 
A s (k) = 2Vl /2 (k)/5 and A T (k) ~ ^ 1/2 

and the running of scalar tilt as at the Hubble radius crossing k 
by: 

dns 



V^ /Z (k)/10, the scalar and tensor spectral indexes ns,T 

aH are defined as usual 



<MA 2 S 
dink 



k=aH 



dlnA 2 ^ 
dink 



a s 



k=aH 



dink 



(27) 



k=aH 



4. The Markov Chain Monte Carlo methodology 

We use the Markov Chain Monte Carlo (MCMC) technique to reconstruct the inflation- 
ary potential and to derive constraints on the inflationary observables in the 4-D inflation 



lr rhe normalization of Vs(k) ensures that A $ coincides precisely w ith the density contrast Sh at Hubble 



radius crossing as defined by Liddle and Lyth (jLiddle fc Lvth 20001 ). The normalization of VtQz) is then 
chosen so that e# = Aj,/Ag. 
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and 5-D single-field i nflation models by using the WMAP 5-year data (IDunkley et al. 2009 
Komatsu et al. 20091 ) complemented with geometric probes from the Type la supernovae 
(SN) distance-redshift relation and the baryon acoustic oscillations (BAO). 
The SN distance-redshift relation has been studied in detail in the re cent unified analysis o f 
the published heterogeneous SN data sets - the Union Compilation08 (IKowalski et al. 20081 ). 
The BAO in the distribution of galaxies are extracted from the Sloa n Digital Sky Survey s 
(SDSS) and Two Degree Field Galaxy Redshidt Survey (2DFGRS) JPercival et al. 2007h . 
The CMB, SN and BAO data (WMAP5+SN+BAO) are combined by multiplying the like- 
lihoods. We decided to use these measurements especially because we are testing models 
deviating from the standard Friedmann expansion. These datasets properly enables us to 
account for any shift of the CMB angular diameter distance and of the expansion rate of the 
Universe. 

For the forecast from PLANCK-like simulated data we use the CMB temperature (T) and 
polarization (P) power spectra of our fiducial cosmological model and t he expected experi- 
mental characteristics of t he Planck frequency channels given in Table 1 (IMandolesi et al. 2009 
Planck Consortia 20051 ) . For each frequ ency channel we consider an homogeneou s detector 
noise with the power spectrum given by (IPerotto et al. 20061 ; iPopa fc Vasile 20071 ): 



I.U 



(# b A a ) 2 exp^ + W ln2 



ce(T,P) 



(28) 



where v is the frequency of the channel, 0& is the FWHM of the beam and A c are the 
corresponding sensitivities per pixel. The global noise of the experiment is obtained as: 



Ew 



(29) 



Our fiducial model is the standard ACMB cosmological model with the physical baryon 
density VL^h 2 = 0.022, the physical dark matter density Q c d m h 2 = 0.11, the ratio of the 
sound horizon distance to the angular diameter distance 6 S = 1.04, the reionization optical 
depth r = 0.085, the scalar spectral index n s = 0.96 and the curvature fluctuations amplitu de 
A 2 S = 2.28 x 10~ 9 at pivot scale k=0.01Mpc~ 1 JDunklev et al. 20091 : komatsu et al. 2009h . 
We present in Fig. 1 the CMB angular power spectra of the the fiducial ACMB cosmological 
model and the temperature and polarization noise power spectra obtained for the Planck 
experimental characteristics presented in Table 1 considering a coverage of the sky of 80%. 
The final evaluation of the systematic effects that could remain in the Planck data after 
data reduction affecting scientific exploitation will come from accurate in-flight analyses and 
extensive Monte Carlo simulations and is out of the scope of this work. On the other hand, 
we include in this study also a degradation of Planck ideal sensitivity possibly introduced 
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by residuals of systematic effects at low multipoles, where the cross-check for systematics 
possible at high multipoles comparing different sky areas is obviously not feasible. The 
most critical source of contamination wil l likely come from the straylight, e.g. the signal 



entering far sidelobes fjSandri et al. 20041 ) at large angular distance from the main beam . 



Two different sources mainl y contribute to this eff ect: the CMB dipole ( jBurigana et al. 20041 ) 



and the Galactic emission (jBurigana et al. 20061 ). The former affects only even multipoles 



but it is larger in amplitudes at the considered frequencies, the latter is smaller in amplitudes, 
but affect all multipoles. Our simple conservative toy model, based on the above studies, 
assumes an increasing of the noise power at low multipoles coming from residuals of the 
angular power spectra estimated for these systematic effects possible generated by a non 
perfect subtraction of them, as in the case in which the properties of Planck optical response 
in the far sidelobes is known only with an accuracy of about 30%. 
The uncertainty added to the instrumental (receiver) noise is clearly visible in Fig. 1. 

We evaluate the likelihood function for 4-D and 5-D inflationary models by using the 



public packages CosmoMC and C AMB ([Lewis & Briddle 200(1 iLewis et al. 1999T ) modified 



to enable us to include the corresponding Hamilton- Jacobi formalism as described in the 
previous section. We perform the analysis in the framework of the flat ACDM standard 
cosmological model. 

For 4-D inflation case, the ACDM standard cosmological model is described by the following 
sets of parameters receiving uniform priors: 

jji2 jjii H'" H' 
Vtjjh , Q c h , S , t , A s , j j 

where H', H" and H'" are the derivatives of the Hubble expansion rate H with respect to the 
scalar field. As notted before by Lesgourgues at al. (2008), because the physical effects in 
the primordial power spectra depend on combinations of Hubble expansion rate derivatives, 
the basis of parameters receiving uniform priors should consists in functions of the above 
combinations or linear combinations of them, ensuring that Markov Chains can converge in 
a reasonable amount of time. 

By analogy, we take for 5-D inflation case the following basis of parameters receiving uniform 
priors: 

Q b h 2 , Q c h 2 , 6 S , r, A 2 S , — , — , , p. 

yl y yl 

where y', y" and y'" are the derivatives with respect to the scalar field of the parameter 
y(4>) defined by Eq.( fT2l . One should note that the 5-D inflation case requires the additional 
parameter /i that controles the hierarchy of 4-D and 5-D Planck mass scales through the 
brane tension A: 

v = — = \hr- 30 

m 4 V 3 
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We run 32 Monte Carlo chains per model and dataset , imposing for each case the Gelman 



& Rubin convergence criterion (IGelman fc Rubin 19921 ) . 



5. The results: analysis and interpretation 

5.1. The 4-D and 5-D inflationary parameter bounds 

The parameter bounds derived from each set of chains are given in Table 2 while Fig. 2 
and Fig. 3 show the results of our fits of 4-D and 5-D inflationary models on WMAP5+SN+BAO 
dataset and PLANCK-like simulated dataset. All parameters are computed at the Hubble 
radius crossing £;*=(). OlMpc -1 . 

From the fit of 4-D inflation model to WMAP5+SN+BAO dataset we obtain bounds on 11$, 
A s , as and R at =0.01 Mpc -1 that translated into bounds at =0.002 Mpc -1 show 
a good agreement wi th bounds on the similar parameters reported by the WMAP team 



(IKomatsu et al. 20091 ). Although our computation does not involve the HSR approximation 
in the computation of perturbations spectra, our constraints on the inflationary parameters 
obtained from the fit of 4-D inflation model to WMAP5+SN+BAO dataset are in general 
in agreement with the similar results obtained by using the HSR formalism to recover the 
i nflationary potential, when imposing constraint s on the number of e-folds during inflation 



(IPeiris fc Easther 2006al ; iPeiris &: Easther 2006bl ). Our results are directly comparable with 



the results presented in Hamann et al. (2008) that uses the same numerical evaluation of 
the spectra to obtain constraints on the inflationary parameter by using a selection of CMB 
data including WMAP complemented by LSS measurements. 

Looking at Fig. 2 and Fig. 3 we see that the amplitude of scalar power spectrum A 2 S ob- 
tained from both datasets is suppressed in the 5-D inflation case, when compared with the 
similar values obtained in 4-D standard inflation. Morever, in 5-D inflation model the joint 
confidence regions of the scalar spectral index n s and tensor to scalar amplitude ratio R are 
anti-correlated, while the amplitude of the tensor power spectrum A\ is increased, when 
compared with the 4-D standard inflation case. This can be attributed to a larger contri- 
bution of the tensor modes to the primordial perturbations in braneworld inflation. In this 
case the normalization of the scalar perturbation is reduced to obtain the correct value of 
R. Likewise, we obtain a strong correlation between the values of R and the tensor spectral 
index tit from the fits of 4-D and 5-D inflation models to both datasets, in good agreement 
with the predictions of consistency relation given in Eq. (1). 

Fig. 4 presents the constraints on HSR parameters derived from the fits of 4-D and 5- 
D inflationary models on WMAP5+SN+BAO dataset and PLANCK-like simulated dataset. 
We show the division of sh-Vh plane into large field, small field and hybrid field classes of 
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inflation models (IKinney 20021 ; iLiddle fc Taylor 20021 ) overlaid with our constraints on their 
joint 68% and 95% confidence intervals. We see that all three classes of inflation models are 
allowed at 2-a level by the fit of 4-D standard inflation to our datasets. 
The joint marginalized distribution of e# and rju obtained from the fit of 5-D inflation 
model shows that the large field and small field classes of inflationary models are allowed 
by WMAP5+SN+BAO and by Planck datasets at 2-a level, while the hybrid class of in- 
flationary models seems to be disfavored by both datasets in the 5-D single-field inflation 
scenario. The parameter values within each class of allowed inflationary models are tightly 
constrained by Planck dataset. 

Of particular interest are the differences between the degeneracy directions in €h-Vh plane 
found from the fit of 4-D inflation model to WMAP5+SN+BAO dataset and Planck dataset 
that arise due to the depe ndence of a. q on £^ . The role of th e £ |f in the dynamics of inflation 
is discussed in details in (jChongchitnan &: Efstathiou 20051 ; lEasther &; Peiris 20061 ) and the 
accuracy of slow-r oll inflation models with signicant running is probed by using Mo nte Carlo 
reconstruction in (lEasther fc Kinney 20031 ; iMakarov 20051 ; iPeiris fc Easther 2006b! ) . 
Looking at Fig. 4 one can see the preference of WMAP5+SN+BAO dataset for large and 
positive £ 2 values in 4-D standard inflation case, that translates into large negative values 
of the running of scalar spectral index a s and a larger degeneracy in ns-us plane, when 
compared to the similar results obtained from the fit to Planck dataset (see Fig. 2). 
The differences between the degeneracy directions obtained in 4-D standard inflation case 
and 5-D inflation arise via the dependence of HSR parameters on the dynamical equations 
driving inflation, which are different in 4-D and 5-D inflationary models. 



5.2. Reconstruction of 4-D and 5-D inflationary potential 



The aim in the reconstruction of the inflationary potential is to take the measurements 
of various inflationary observables corresponding to a particular wavenumber k and to use 
them to obtain the inflationary potential V(4>) and its derivatives at the scalar field value 
0* when the scale k crosses the Hubble radius k* during inflation. 

In the general case of the single-field braneworld inflation, the slope and the curvature of 
the 5-D inflationary potential as function of inflation ary observables R and n s and on the 
combination V/ A are given by (ILiddle fc Taylor 20021 ) : 



V 

Yl 



16ttR 



1 + 



1 + V/2X 
L G(V/X) J 
V 



2A 



6R 



l + V/X 
G 2 (V\) 



+ (n s - 1) 



(31) 
(32) 
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where: 



G(x) 



yl + x 2 — x 2 sinh - 



x 



-1/2 

(33) 



In the high-energy limit (V » A) the function G 2 (V/X) — > 3V/2X. In the low-energy limit 
(A >> V) G 2 {y/\) — > 1 and the scalar and tensor perturbation spectra of 4-D standard 
inflation are recovered. 

We compute the magnitude, the slope and the curvature of the inflationary potential 
from the fits of 4-D and 5-D inflation models to our datasets by using Eq.(lO) and Eq.(18) 
respectively. In Fig. 5 we show the allowed regions of the recovered magnitude of the in- 
flationary potential, its slope and curvature from the fit of 4-D and 5-D inflation models 
to WMAP5+SN+BAO dataset and PLANCK-like simulated dataset. We also show the ID 
marginal distribution of the recovered inflationary potential, its slope and curvature as a 
function of V/ A, obtained from the fit of 5-D inflation model to the same datasets. 
Fig. 5 explicitly demonstrates the effect of the braneworld reconstruction of the inflationary 
potential. As Vj A is increased the magnitude and the curvature of the inflationary potential 
are decreased while its slope steepens. Also, the magnitude, the slope and the curvature of 
the inflationary potential are increased in both 4-D and 5-D inflationary scenatios when R 
increases. 

The results from the fit of 5-D inflation model to WMAP5+SN+BAO dataset and PLANCK- 
like simulated dataset show that the magnitude and the slope of the inflationary potential 
are anti-correlated with the scalar spectral index, n s . 

The mean values of the magnitude, slope and curvature of the inflationary potentials to- 
gether with their 95% upper and lower intervals are given in Table 2. 

In Fig. 6 we show the dependence of the reconstructed regions of 4-D and 5-D inflationary 
potentials allowed by the same datasets (at 65% CL) in an observational inflationary window 
corresponding to AN = 11 e- folds, as functions of the scalar field. 



5.3. The 4-D and 5-D single- field inflation consistency relations 



Th ere is an infinite h i erarchy of consistency equations of the single-field standard in 



flation (ILidsev et al. 19971 : iSong k Knox 20031 : IChung et al. 20031 : iChung k Romano 2006 
Cortes k Liddle 2006 I ). To the leading order in the slow-roll approximation, the consistency 



relation of the standard scenario given in Eq.(l) is degenerate. To next-tq-leading order, this 
consistency relation receives corrections of the form (ICopeland et al. 19941 : ILidsev et al. 19971 ) 



Tlx 



A 2 
A 2 



A 2 
A 2 



n s ) 



(34) 
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that do not depend on the spectral index of the tensor perturbations. 

As the inflationary observables 71$, and R are evaluated at the epoch of horizon-crossing 
quantified by the number of e-folds iV before the end of the inflation at which our present 
Hubble scale equalled the Hubble scale during inflation, the uncertainties in the determina- 
tion of N translates to theoretical errors in the determination of the inflationary observables. 
Assuming that the ratio of the entropy per comoving interval today to that after reheating is 
negligible, the main uncertainty in the determination of N is caused by our ignorance in the 
determination of the reheating temperature after inflation leading to an error of AiV ~ 14 
fcinnev fc Riotto 20061 : Adshead fc Easther 2008^ . 

In order to test the observational signature that standard and braneworld inflationary sce- 
narios may produce, we use the estimates of the inflationary parameters obtained from the 
fits to WMAP5+SN+BAO and PLANCK-like simulated datasets to compare the experimen- 
tal difference between tensor spectral indexes, n^P — n^P , to the theoretical error in the 
tensor spectral index computed by using the consistency relation fl34"|) . 
To the lowest order in slow-roll parameters, the uncert ainties AR and Ans in terms of the 
uncertainty in the number of e-folds AiV are given by (IKinney 20021 ; iKinney et al. 20041 ): 



AR 
AN 



R 



(ns ~ 1) 



R 



An s 
AN 



v. 



-r 2 + ie 

32 ^ 



(35) 
(36) 



The theoretical uncertainty on the tensor spectral index in the standard 4-D inflation can 
be straightforward obtained from Eq. (l3~4|) by using Eqs.(35) and (36): 



An c T 
AiV 
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r 



1 _ ns_ _ A 
2 Al 



AR 
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+ 



1 A\ An s 
8 A| AiV 



(37) 



The estimate of An^/AN should be compared to An? = n. In Table 3 we present the 
mean values of the lowest order estimates of the theoretical errors AR/R, Ans/ns and An^ 
from the fit of 4-D inflation model to WMAP5+SN+BAO and Planck datasets obtained 
by assuming AiV=14 e-folds and the mean values of the difference between the experimen- 
tal values of the tensor spectral indexes An? = n^P — nip, while in Fig. 7 we show their 
ID marginalized likelihood probability distributions. We also show in the same figure the 
ID marginalized likelihood probability distribution of the lowest order estimates of the the- 
oretical errors An^,, obtained by assuming AiV=14, compared with the ID marginalized 

and their 2D joint 



n, 



ID 

r 



5D 



likelihood probability distribution of the difference An^ 
allowed bounds (at 68% and 95% CL). 

The analysis of the results presented in Fig. 7 and Table 3 shows that the Aht parameter 
space obtained from the fits of 4-D and 5-D inflation models to WMAP5+SN+BAO dataset 
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is dominated by the theoretical error An^: the confidence interval corresponding to An? is 
smaller by a factor of 1.2 than that corresponding to Any. The same parameter space is 
better constrained by the Planck dataset: in this case the confidence interval correspond- 
ing to An^ is three times smaller than that corresponding to Aht- 

We conclude that the detection of tensor perturbations and the theoretical uncertainties in 
the inflationary observable represent a significant challenge for the future Planck CMB 
measurements: distinguishing between the observational signatures of the standard and 
braneworld single-field inflation scenarios. 



6. Conclusions 

One of the most anticipated results of forthcoming Planck high precision CMB mea- 
surements is probing the physics of inflation and in particular, the reconstruction of the 
inflation potential. On the other hand, the possibility that our four-dimensional Universe 
could lie on a brane embedded in a higher dimensional space has important consequences 
for the early universe and in particular for the cosmological inflation. 

In this paper we make a more general determination of the inflationary observables in 
the 4-D and 5-D single-field inflationary scenarios by exact reconstruction of the dynamics 
of the inflation potential during the observable inflation, with a minimal number of assump- 
tions. 

Making use of the general forma l ism for 5-D single-field in flation developed by Hawkins and 



Lisdey (jHawkins &: Lidsey 200ll ; lHawkins Lidsey 20031 ) valid in all regimes, having many 



of the properties of the Hamilton- J acobi formalism developed for the 4-D standard inflation, 
we compute the scale dependence of the amplitudes of the scalar and tensor perturbations 
by integrating the exact mode equation. Our computation does not assume the slow-roll 
approximation and is valid in all regimes if the field is monotonically rolling down its po- 
tential. The solutions in 4-D and 5-D inflation scenarios differ through the dynamics of 
the background scalar field and the number of e-folds assumed to be compatible with the 
observational window of inflation. 

We address higher-order effects in the standard and braneworld single-field inflation scenar- 
ios by fitting the Hubbble expansion rate H(<p) and subsequently the inflationary potential 
V ((/)), directly to WMAP5+SN+BAO and Planck-like simulated datasets. 
One should note that our results refer to the initial scalar and tensor perturbation spectra 
and not to the braneworld eff ects on the subsequent evolution of th e perturbations that is 



likely to be model dependent (ILeong et al. 20021 ; iRhodes et al. 20031 ) 



Assuming that the ratio of the entropy per comoving interval today to that after reheating 
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is negligible, we analyze the implications of the theoretical uncertainty in the determination 
of the reheating temperature after inflation on the observable predictions of inflation. 
We find that the detection of tensor perturbations and the theoretical uncertainties in the in- 
flationary observables represent a significant challenge for the future Planck CMB measure- 
ments: distinguishing between the observational signatures of the standard and braneworld 
single-field inflation scenarios. 

We acknowledge the use of the GRID computing system facility at the Institute for 
Space Sciences Bucharest and to the staff working there. 

L.P and A.C. are partially supported by ES A/PECS Contract C98051 and CNCSIS Contract 
539/2009. 



Table 1. The expected experimental characteristics for the Planck fre quency channels 
considered in the paper (IMandolesi et al. 20091 ; IPlanck Consortia 20051 ) . At and A p are 
the sensitivities per pixel for temperature and polarization maps. 



Frequency (u) FWHM A T A P 
(GHz) (arc-minutes) (fiK) (;uK) 



70 

100 

143 



13 
9.5 
7.1 



23.48 33.21 
6.8 10.9 
6.0 11.4 
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Table 2. The mean values and 95% CL lower and upper intervals of the derived 
parameters obtained from the fit of 4-D and 5-D inflation models to WMAP5+SN+BAO 
dataset and PANCK-like simulated dataset. All parameters are computed at the Hubble 

radius crossing A;*=0.01 Mpc -1 . 



WMAP5+SN+BAO Planck 
Parameter 4-D Inflation 5-D Inflation 4-D Inflation 5-D Inflation 
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Fig. 1. — The CMB angular power spectra of the fiducial ACMB cosmological model and the 
temperature (dashed red line) and polarization (dash-dotted blue line) noise power spectra 
obtained for the Planck experimental characteristics presented in Table 1, considering a 
coverage of the sky of 80%. 



Table 3. The mean values and 95% CL lower and upper intervals of the lowest order 
estimates of the theoretical errors AR/R, Aris/ris and An^ obtained from the fit of 4-D 
inflation model to WMAP5+SN+BAO and Planck datasets by assuming AA/=14, and 
the mean values and 95% CL lower and upper intervals of An? = n^P — rijP, obtained from 
the fits of 4-D and 5-D inflation models to the same datasets. 
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Fig. 2.— We show the results of the fits to the WMAP5+SN+BAO dataset of 4-D inflation 
model in red and 5-D in inflation model in blue. The results of the fits to the PLANCK-like 
simulated dataset of 4-D inflation model are in magenta and of 5-D in inflation model are 
in cyan. The top plot in each column shows the probability distribution of different scalar 
inflationary observables while the other plots show their joint 68% and 95% confidence 
intervals. All parameters are computed at the Hubble radius crossing fc*=0.01Mpc _1 . 
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Fig. 3.— We show the results of the fits to the WMAP5+SN+BAO dataset of 4-D inflation 
model in red and of 5-D in inflation model in blue. The results of the fits to the PLANCK-like 
simulated dataset of 4-D inflation model are in magenta and of 5-D in inflation model are in 
cyan. The top plot in each column shows the probability distribution of different tensorial 
inflationary observables while the other plots show their joint 68% and 95% confidence 
intervals. All parameters are computed at the Hubble radius crossing fc^=0.01Mpc _1 . 
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Fig. 4. — We show the bounds on the HSR parameters derived from the fits to 
WMAP5+SN+BAO dataset of 4-D inflation model in red and of 5-D inflation model in 
blue. The results of the fits to the PLANCK-like simulated dataset of 4-D inflation model are 
in magenta and of 5-D in inflation model are in cyan. The top plot in each column shows the 
probability distribution of different HSR parameters while the other plots show their joint 
68% and 95% confidence intervals. All parameters are computed at the Hubble radius cross- 
ing /c*=0.01Mpc -1 . We show the division of e# — t]h plane into large field (— e# < Vh < 
small field (t]h < — e#) and hybrid field (0 < tu < Vh) classes of inflation models. 
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VA. R n s 

Fig. 5. — We presesent the recovered regions (at 68% and 95% CL) of the magnitude, slope 
and curvature of the inflationary potential from the fits to WMAP5+SN+BAO dataset of 
4-D inflation model in red and of 5-D inflation model in blue. The similar results from the 
fits to PLANCK-like simulated dataset of 4-D inflation model are in magenta and of 5-D in 
inflation model are in cyan. We also show the ID marginal distribution of the recovered 
magnitude of the inflationary potential, its slope and curvature (at 95% CL) as function 
of V/X obtained from the fit of 5-D inflation model to WMAP5+SN+BAO dataset (black 
lines) and PLANCK-like simulated dataset (yellow lines). All parameters are computed at 
the Hubble radius crossing /c^O.OlMpc -1 . 
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Fig. 6. — We present the regions of 4-D and 5-D inflationary potentials allowed at 65% CL by 
WMAP5+SN+BAO and PLANCK-like simulated datasets in the same observational inflation 
window corresponding to AN = 11 e-folds. For each case, the magnitude of inflationary 
potential is normalized to V*, the value of the inflationary potential at Hubble radius crossing 
fe,=0.01Mpc _1 . 
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Fig. 7. — Panels a) and b): ID marginalized likelihood probability distributions of the the- 
oretical errors AR/R and Ans/ns obtained from the fits of 4-D standard inflation model to 
WMAP5+SN+BAO dataset (continuous line) and Plaanck dataset (dashed line) assuming 
AiV=14. Panels c) and d): 2D joint marginalized probability distributions (at 68% and 95% 
CL) obtained from the fits to WMAP4+SN+BAO dataset. The red filled contours corre- 
sponds to Aht = ArijP — ArVjP and the white filled contours corresponds to the theoretical 
error An^ obtained by assuming AiV=14. Panel e): ID marginalized likelihood probability 
distributions of An9p (black continuous line) and An? = n^P — n^P (red continuous line) 
obtained from the fits to WMAP5+SN+BAO dataset. Panels e) and g): 2D joint marginal- 
ized probability distributions (at 68% and 95% CL) obtained from the fits to PLANCK-like 
simulated dataset. The magenta filled contours corresponds to Any = An^P — An^P and the 
white filled contours corresponds to the theoretical error An^ obtained by assuming AiV=14. 
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